Noise and Signal scaling factors in Digital Holography in weak illumination: 

relationship with Shot Noise. 
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We have performed off axis heterodyne holography with very weak illumina- 
tion by recording holograms of the object with and without object illumination 
in the same acquisition run. We have experimentally studied , how the recon- 
structed image signal (with illumination) and noise background (without) scale 
with the holographic acquisition and reconstruction parameters that are the 
number of frames, and the number of pixels of the reconstruction spatial fil- 
ter. The first parameter is related to the frequency bandwidth of detection 
in time, the second one to the bandwidth in space. The signal to background 
ratio varies roughly like the inverse of the bandwidth in time and space. We 
have also compared the noise background with the theoretical shot noise back- 
ground calculated by Monte Carlo simulation. The experimental and Monte 
Carlo noise background agree very well together. 
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1. Introduction 



Demonstrated by Gabor^ in the early 50's, the purpose 
of holography is to record, on a 2D detector, the phase 
and the amplitude of the light shining from an object un- 
der coherent illumination. Since a thin film does not pro- 
vide a direct access to the recorded data, the holographic 
film has been replaced by 2D electronic detector in dig- 
ital holography^, whose main advantage is to perform 
the data acquisition and the holographic reconstruction 
numerically^^^. Off- axis holography^ is the oldest and the 
simplest configuration adapted to digital holography^^^^^. 
In off-axis digital holography, as well as in photographic 
plate holography, the reference or local oscillator (LO) 
beam is angularly tilted with respect to the object ob- 
servation axis. It is then possible to record, with a single 
hologram, the two quadratures of the object's complex 
field. However, the object field of view is reduced, since 
one must avoid the overlapping of the image with the 
conjugate image alias. Phase-shifting digital holography, 
which has been later introduced^, records several images 
with a different phase for the LO beam. It is then possi- 
ble to obtain the two quadratures of the field in an in-line 
configuration even though the conjugate image alias and 
the true image overlap, because aliases can be removed 
by taking image differences. 

With the development of CCD camera techniques. 
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digital holography became a fast-growing research field 
that has drawn increasing attention^^^^. Off- axis 
holography has been applied recently to particle^, 
polarizationi^, phase contrast synthetic aperturei^, 
low-coherence^^^^^ and microscopic^^^^ imaging. Phase- 
shifting holography has been applied to 3D^^^^^, color—, 
synthetic aperture^^, low-coherence^^, surface shaped! 
and microscopic^^'^^ imaging. 

We have developed an alternative phase-shifting dig- 
ital holography technique that uses a frequency shift of 
the reference beam to continuously shift the phase of the 
recorded interference pattern^^. 

More generally, our setup can be viewed as a multi- 
pixel heterodyne detector that is able of recording the 
complex field E scattered by an object in all pixels of the 
CCD camera in parallel. Because the holographic signal 
results in the interference of the object signal complex 
field E with a reference (LO) complex field Elo^ whose 
amplitude can be much larger (i.e. \Elo\ ^ |^|), the 
holographic detection benefits of "heterodyne gain" (i.e. 
\E E^q\ ^ |£^P), and is thus well suited to detect weak 
signal fields E. 

By performing phase shifting heterodyne holography 
in an off-axis configuration, the holographic signal can be 
frequency shifted both in time and space. It is then pos- 
sible perform a very efficient filtering of the holographic 
signal in space and time. This double filtering combined 
with heterodyne gain open the way to holography of weak 
object fields E with ultimate sensitivity^^. This type 
of holography has been used (i) to detect and analyse 
the laser Doppler spectrum of the extremely weak opti- 
cal signal that is transmitted and scattered by a human 
breast in vivo (< 1/10 photoelectron per pixel) (ii) to 
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image, localize and analyse the motion of metallic nano 
particules^^'^^ , and (iii) to analyse vibration of very small 
amplitude^"— by detecting the holographic signal on the 
vibration sideband^^. 

Since at low signal the prevalent and the limiting noise 
is shot noised, it is important to study how shot noise 
limits the image quality in digital holography. Charriere 
et al., have studied, with simulated and experimen- 
tal data, how shot noise influences the quality of the 
phase images, in particular in the case of short expo- 
sure time^^^^^. On the other hand, Verpillat et al. have 
calculated, in the weak illuminated case, the shot noise 
limits for amplitude images^^, and have shown that the 
shot noise background can be easily calculated by Monte 
Carlo simulations. 

In this paper, we have considered amplitude images by 
performing off axis heterodyne holography in very weak 
illumination conditions. We have modified our experi- 
mental setup to record holograms with and without illu- 
mination of the object (but with reference beam) in the 
same acquisition run. We have measured the noise back- 
ground without illumination and compared it with the 
signal obtained with illumination. We have studied ex- 
perimentally, how the signal and the noise scales, on am- 
plitude reconstructed images, with the two mains param- 
eters of the holographic acquisition and reconstruction 
that are the number of frames, and the number of pixels 
of the reconstruction spatial filter—. The first parameter 
is related to the detection frequency bandwidth in time, 
the second one to the bandwidth in space. The expected 
scale factor values are retrieved: the signal to noise back- 
ground ratio varies roughly like the inverse of the band- 
width in time and space. Finally, we have compared the 
noise background observed in experiments, with the shot 
noise background calculated by the Verpillat et al. Monte 
Carlo method. The experimental and Monte Carlo noise 
background agree very well together. This demonstrates 
that our weak illumination holography is shot noise lim- 
ited. 

2. Experimental Setup and Numerical Recon- 
struction 

The experimental setup is shown in Fig. [TJa). The light 
source is a A = 532 nm, 50 mW laser (SLIM-532: Oxxius 
SA). The laser beam (complex field El, frequency cjl) is 
split into an illumination beam {Ej, uji) and a reference 
local oscillator (LO) beam {Elo^ ^lo)- The LO beam 
is enlarged by a beam expander BE, and is tilted by an 
angle 0^1° with respect to the observation axis of the 
CCD camera. 

The illumination beam is enlarged by a short focal 
length (5 mm) objective O, located at a distance D'—D 
63 cm in front of an United States Air Force (USAF) tar- 
get. The light transmitted by this object is superimposed 
on the LO beam. The interference pattern is recorded by 
the CCD camera (PCO Inc. Pixelfly digital camera: 12 
bit, frame rate ujccd=^'^-^ Hz, exposure time T = 50 
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Fig. 1. (a) Digital holography setup. L: diode laser; 
BSl, BS2: beam splitters; AOMl and A0M2: acousto- 
optic modulators; O: objective; BE: beam expander; M: 
mirror; A: light attenuator. USAF: transmission USAF 
target to image. CCD : CCD camera, (b) Illumination 
beam amplitude (solid line) and phase (dashed line) as a 
function of the frame index. 



ms, with 1280 x 1024 pixels of 6.7 x 6.7 /im) located 
dit D = 32 A cm downstream. The signal of the camera 
measured in Digital Count (DC) units (0 to 4095) is mul- 
tiplied by the camera gain {G = 4.41 e/DC as calibrated 
by PCO Inc.) to be converted in photoelectron (e) units. 
A set of optical attenuators A (gray neutral filters) is 
used to reduce the intensity of the illumination beam. 

Two acousto-optic modulators (Crystal Technology: 
^AOMi,2 — 80 MHz) are used to adjust ujj and ujlo- 



(1) 



such way the phase of the EE^q term is shifted by 7r/4 
from one CCD frame to the next (8 phases detection). 
We have thus: 

^LO — UJi = UJAOM2 — ^AOMl = ^CCd/^ (2) 

In addition, the ujaomi — 80 MHz signal that drives 
AOMl is modulated in amplitude at ujccd/2 through 
the switch SW, so that the illumination beam is turned 
on and off from one CCD frame to the next. The am- 
plitude and phase (relative to the LO beam) of the illu- 
mination beam are displayed in Fig(T] (b). Even frames 
are recorded with the illumination beam on the object, 
whereas odd frames are recorded without. We have then 
recorded sequences of 64 successive CCD frames (/q to 
/ea) that are analyzed to reconstruct the intensity images 
of the USAF target. 

Compared with the heterodyne holography setup of 
ref^^, the Fig(T]setup exhibits two main differences. First, 
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the AOMs acts on both the niumination and local oscil- 
lator arms (and not on the local oscillator arm only). It 
is thus possible to put the illumination attenuator A up- 
stream with respect to AOMl. The parasitic light (of fre- 
quency ujl) that is diffused by the optics upstream with 
respect to the attenuator A is then frequency shifted by 
80 MHz with respect to the detected frequency ( ulo)- 
This parasitic light is thus not detected, and the illumi- 
nation leaks are reduced. 

Second, the electronic switch turns on and off the illu- 
mination in such a way that the even frames are recorded 
with illumination of the sample, and the odd frames with- 
out (but the LO beam remains for both even and odd 
frames). In the following "with illumination" means thus 
with illumination and LO beams on, while "without il- 
lumination" means with illumination beam off and LO 
beam on. The frequency shift of the LO is here ojccd/^ 
(and not uccd/^) to record standard 4 phases hologram 
with and without object illumination: the phase shift 
between successive even (or odd) frames is therefore 7r/2. 

We must notice that the combination of the electronic 
switch with the 8 phases detection is not needed for mak- 
ing 4 phase phase shifting holography. It has been used 
here to be able to extract from the same sequence of data 
different types of holograms. With the even frames (or 
odd frames), one can perform 4 phase holography with 
or without illumination of the sample. This can be done 
with 4, 8, 16 or 32 frames. One can also consider only 
the 0^^ and 4^^ (or 1^^ and 5^^) frames to make 2 phase 
holography with or without object illumination. By us- 
ing the 0^^ frame and making the difference with the 1^^ 
frame one can also perform single shot off-axis hologra- 
phy with suppression of the DC term (as proposed by 
Kreis et al.^^). In the following we have analyzed the re- 
constructed images corresponding to these different types 
of hologram to better understand how the signal and the 
noise depend on the number of frames. 

Many numerical methods can be used to reconstruct 
the image of the object^^^^^. The most common methods 
are the Schnars et al. methods^ that involves a single 
discrete Fast Fourier Transform (FFT), and the angular 
spectrum method that involves two FFT^^^^^. In the last 
case, it is possible to use the Cuche et al. method^^ to 
perform a spatial filtering of the pertinent holographic 
signal. 

We have consider here the angular spectrum method 
with spatial filtering, and we have studied how the holo- 
graphic signal and noise scale with the number of frames, 
and with the area of the spatial filter. 

To illustrate the reconstruction (see Fig. [2j), we have 
considered a phase shifting complex hologram H with 4 
phases: H = {Iq — I2) + j(/i — /3), where /q, h ....I3 are 
the CCD signals recorded on 4 successive frames and = 
— 1. The intensity of the CCD plane complex hologram 
|i^(x,?/,0)p (where x and y are the coordinate of the 
CCD pixels) is displayed on Figl2] (a). The reciprocal 
space hologram H{kx^ky^O) (where kx and ky are the 
coordinates of the reciprocal space) is calculated from 




Fig. 2. (a) Intensity signal \H{x,y^z = 0)p de- 
tected by the CCD camera. (b,c) k-space intensity signal 
\H{kx^ky^ z = 0)p without (b) and k-space filtering (c) 
with two filter sizes (400 x 400. (d) Reconstructed im- 
ages \H{x,y,z = D)\'^ of the USAF target. Display is in 
arbitrary log scale. 



the CCD plane hologram H{x^ z = 0) by a Fast Fourier 
Transform (FFT): 



H{kx,ky,0) = FFT [H{x,y,0)] 



(3) 



The reciprocal space hologram H{kx,ky^O) is displayed 
on Fig. [2] (b). The USAF target information lies in 
the bright zone in the left hand side of Fig. [2] (b) that 
corresponds to order +1 image. Because of the 4 phase 
detection, the zero order and order -1 images are lower in 
intensity. There are nevertheless still visible because of 
logarithmic scale display. To select the order +1 signal, 
we have here performed spatial filtering with a 400 x 400 
pixels filter. A 400 x 400 pixels square region centered on 
order +1 is thus cropped, and translated to the centre of 
a zero-filled 1024 x 1024 calculation grid (zero padding) 
as shown on Figl2fc). The cropping operation performs 
a spatial filtering in the Fourier spaced, while the trans- 
lation compensates for the off- axis geometry. 

In the reciprocal space, the hologram H is then mul- 
tiplied by the quadratic phase factor that describes the 
optical field propagation from the CCD plane z = to 
the USAF target plane z = D 

H{kx^ky,z)=H{kx^ky,0) e^i^"^^>/^o^ (4) 



where ko = y fc^ + + = 27r/A. The reconstructed 
image H{x, y, D) is then obtained by inverse FFT: 



H{x, y, z) = FFT-i \H{k,,ky, z) 



(5) 
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H{x,y,D) is displayed on Figl2fd). Because the recon- 
struction involves 2 FFT, the pixel size is the same than 
in the CCD hologram i.e. 6.7 x 6.7 /im. 

3. Experimental results 

We have recorded series of 64 frames for different illu- 
mination levels of the USAF target, and we have recon- 
structed holographic images of the target by varying (i) 
the parity of the frames, (ii) the number of frames, and 
(iii) the area of the spatial filter. 

To analyze the results, we have quantified the noise 
in different ways. First, intensity images have been dis- 
played in order to get a qualitative idea of the noise. To 
make a more qualitative analysis, we have performed a 
cut along a line crossing the black bars of the USAF tar- 
get (like line AA' of Figl2] (c)), and we have plotted the 
intensity along this cut. We have averaged the pixel in- 
tensity over the whole reconstructed image, and we have 
plotted this average intensity as a function of the recon- 
struction parameter (parity, number of frames, and area 
of the spatial filter). We have finally compared the noise 
background with the one resulting from shot noise. 

A. Images and cuts 

We have first studied how the signal, and the noise de- 
pends on the number of frames recorded with and with- 
out signal. From the sequence of CCD frames In (with 
n = to 63), we have built n frame holograms Hn with 
illumination (even frames) defined by: 

H2 = lo- h (6) 

k=K-l 

H^K = ^ {hk — hk+4) + j(^8fe+2 — hk+e) 
k=0 

with = 2,4,8. In this equation II2 a 2 phase phase- 
shifted hologram and a 4 phase holograms. i^8, ^16 
and IIs2 are 4 phase holograms recorded with 8, 16, and 
32 frames. 

Similarly to the holograms constructed with n even 
frames, we have constructed the holograms H'^ with n 
odd frames (i.e. without object illumination) by replac- 
ing in Eq. ([6]) the even frames /2n by the odd frames 
hn-\-i' We have also consider the single illuminated 
frame hologram Hi^i obtained by subtracting the refer- 
ence frame Ii recorded without illumination of the object 
to the frame /q recorded with illumination (the two in- 
dexes 1 and 1 correspond here to the number of even and 
odd frames)^^: 

Hi,i = lo-h (7) 

Figure [3] shows reconstructed images of the USAF tar- 
get illuminated at very low level (estimated to be 0.16 
photon electron per pixel and per frame) with a spatial 
filter of 282 x 282 pixels. Figure [3] (a) is obtained from 
Hi I (i.e. with one frame with and one frame without 




Fig. 3. Intensity image (i.e. \H{x,y,D)\^) of the USAF 
target reconstructed from Hi i (a), H2 (b), II4 (c), Hg 
(d). Reconstruction with spatial filter 282 x 282 pixels. 



object illumination). The Signal to Noise Ratio (SNR) 
of the holographic image is low, and the black bars of 
the target are barely visible (here, and in the remaining 
of this article, the term SNR must be understood, in a 
very qualitative way: SNR qualify the visual quality of 
either an image or a curve). Figure [3] (b,c) and (d) show 
the reconstructed images obtained from i^2, II4 and Hs 
i.e. with 2, 4 and 8 illuminated frames. As expected the 
SNR is becoming greater and greater. Figure H] is similar 
to FiglSl It shows the reconstructed images at a slightly 
higher illumination level (i.e. 1.8 photon electrons per 
pixel and per frame). As expected, SNR is higher. The 
bars are visible on all images i.e. on Fig. |l](a) to (d). As 
expected the SNR is becoming greater and greater from 
FigH] (a) to (d) i.e. from 1 to 8 frames with illumination 
or from Hi i to Hg. 

In order to perform a more quantitative analysis of 
these results, we have plotted, in logarithmic scale, the 
intensity of the reconstructed image signal (i.e. |i^P) 
along the vertical line AA'. Out of the bars, we expect 
to get a higher signal with some variation, because 
the laser illumination is not uniform, and because illu- 
mination level and SNR are low. Within the black bars, 
the signal must go down to zero, and one expects 
to see grooves corresponding to the black bars. In or- 
der to decreases the noise of the curves, we have aver- 
aged over 10 pixels along the x direction. Since 
the bars are horizontal and since the width of the bar 
is much larger than 10 pixels, we expect to keep roughly 
the same contrast for the grooves than without averaging 
(but with less noise). In order to visualize the intrinsic 
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Fig. 4. Same as Figl3] with higher inumination: inten- 
sity image (i.e. \H{x^y, D)\'^) of the USAF target recon- 
structed from Hi^i (a), H2 (b), H4 (c), Hg (d). Recon- 
struction with spatial filter 282 x 282 pixels. Vertical axis 
units are e^ per pixel. 



noise of the holographic detection, we also have plotted 
in light grey the cuts obtained with the same number of 
frames recorded without object illumination (i.e. with 
odd frames). Figure [5] shows the plots obtained from the 
Figl3]data, i.e. from the frames recorded at very low level 
of illumination. Figure[6]is similar to Figl5]and shows the 
plots obtained from the FigHJdata, i.e. from the frames 
recorded at slightly higher illumination. 

Plot of Fig. [5] (a) is obtained with Hi^i i.e. with 
one frame with illumination, since Hi i is obtained from 
the difference of two frames (/q and /i). We have also 
plotted the cut obtained with H2 i.e. with two frames 
without object illumination (i.e. with the same number 
of frames that for i^i,i). The holographic signal is ex- 
tremely low, and the vanishing grooves corresponding to 
the black bars are not visible. Nevertheless, the curve ob- 
tained from Hi^i is slightly above the one obtained from 
i^2- Plot of Fig. [5] (b) is obtained with H2 i.e. with two 
frames with illumination. To visualize the noise, the plot 
obtained with H2 is also drawn. Since the SNR of the 
image is higher, the grooves from the black bars are now 
visible. Out of the grooves the curve obtained with illu- 
mination is "above" the curve obtained without. Within 
the grooves, the curves obtained with and without object 
illumination are roughly at the same level. Plots that are 
displayed on Fig. [5](c) and (d) are similar to those of Fig. 
[5](b). Out of the grooves the curves obtained with illu- 
mination are "above" the curves obtained without, and 
since the number of frames increases \H\'^ becomes higher 
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Fig. 5. Intensity (i.e. \H{x^y^ D)\'^) along the AA' white 
line of FiglSl The black curves correspond Hi^i (a), H2 
(b), H4 (c) and Hg (d). The light grey curves in back- 
ground correspond to H2 (a), H2 (b), H'^ (c) and (d). 
Reconstruction is made with a spatial filter of 282 x 282 
pixels. Vertical axis units are e^ per pixel. 
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Fig. 6. Same as Fig E] with higher illumination: intensity 
(i.e. |i^(x,?/,i:>)p) along the AA' white line of FigH The 
black curves correspond Hi^i (a), H2 (b), H4 (c) and Hs 
(d). The light grey curves in background correspond to 
H2 (a), H2 (b), and (d). Reconstruction is made 
with a spatial filter of 282 x 282 pixels. 



and higher. Within the groove the curves obtained with 
an without object illumination are at roughly the same 
level. A good estimate of the holographic detection noise 
can be thus obtained by comparing the curve obtained 
with illumination, with the curve obtained in the same 
conditions (same number of images, same spatial filter) 
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without object illumination. 

Plots of Fig. [6] are similar to the ones of Fig. [5l but 
with higher illumination. As expected \H\^ is higher out 
of the grooves. Within the grooves, \H\^ decreases, and 
the relative contrast of the groove becomes higher from 
Fig. [6ja) to Fig. [6l^d). Nevertheless, since the holo- 
graphic signal becomes higher and higher, do 
not reach, within the groove, the noise floor that corre- 
spond to the data recorded without object illumination. 



B. Average intensity \H\^ 

In order to determine how the signal and the noise scale 
with various experimental and reconstruction parameters 
(parity of the frames, number of frames, and area of the 
spatial filter), we have averaged the pixel intensity over 
the whole reconstructed image, and we have plotted, for 
the even and odd frames, the average intensity \H\^ as a 
function of either the number of frames, or the area of the 
spatial filter. The average intensity obtained with illu- 
mination is then mainly related to the holographic signal 
from the USAF target, while the intensity obtained with- 
out object illumination is related to the detection noise. 
To simplify the analysis, we have considered the case of 
the higher illumination level of FigH] and Figl6fa). In 
that case, the grooves can be seen with only one illumi- 
nated frame (i.e. for as seen on FiglHl This means 
that signal from the USAF target overcomes the noise in 
the worst case. 

Figure [3 shows the averaged (over the 1024 x 1024 pix- 
els of the reconstructed images) intensity \H\^ . This 
intensity is plotted for both even and odd frames as a 
function of the number A/'^ of frames recorded with illu- 
mination in Fig. [Tl^a) or as a function of the area of the 
spatial filter in Fig. [Tl^b). 

Consider first Fig. [7] (a). For even frames (i.e. with il- 
lumination) the average intensity \H\^ increases quadrat- 
ically with the number of frames (with a Log Log slope 
close to 2). This is expected, since we perform a coher- 
ent acquisition of the holographic data. The amplitude 
of the holographic signal (i.e. H) increases linearly with 
acquisition time, i.e. linearly with A'5. The intensity |i7 p 
increases thus like N'^. We must notice that the slope is 
slightly lower that 2. This may be due to loss of coher- 
ence of the holographic detection for long time (i.e. large 
N})) and to noise for short time (i.e. small A'fo), since the 
noise contribution becomes larger when the detected sig- 
nal goes down. Note that the noise contribution strongly 
depends on spatial filter area. As a consequence, the av- 
erage energy varies quite a bit for the first point of the 
curves (A^^ = 1 i.e. Hi i ) as seen on Fig. [71(a). 

With odd frames (i.e. without illumination), the aver- 
age intensity \H\^ increases linearly with the number of 
frames A'5 with slope equal to 1, as expected. Since the 
phase of the noise changes from one frame to the next, 
we perform an incoherent acquisition of the noise. The 
amplitude H of the holographic noise increases thus like 
the square root of the acquisition time, i.e. like V^A^. 
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Fig. 7. Average intensity \H\^: (a) as a function of the 
number A'^ of frames recorded with illumination: A'5 = 1 
for or 1, A^b = 2 for H2 or i^^,..., and AT^ = 32 for 
i^32 or i^32, and (b) as a function of the area of the spatial 
filter (in pixel units). Even frames curves are black, odd 
frames curves are light grey. 



The intensity \H\^ increases thus like A'^. 

Figure [7] (b) analyses the dependence of the average 
intensity with the spatial filter area. For ev en fr ames 
(i.e. with illumination) the average intensity \H\^ does 
not strongly depend on the spatial filter area: the black 
curves are flat. This is expected since most of the holo- 
graphic signal energy is within the low spatial frequency 
modes. For odd frames contrarily (i.e. without illumi- 
nation), the average intensity \H\^ varies linearly with 
the spatial filter area: the light grey curves of Fig. [71 (b) 
have a slope equal to 1. This is expected since we de- 
tect noise. The noise is incoherent and homogeneous in 
the spatial frequency region that is selected by the spa- 
tial filter. This point is illustrated by Fig. [8] that shows 
a reciprocal space intensity image H{kx^ky^O) obtained 
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Fig. 8. Reciprocal space intensity image H{kx^ ky^O) ob- 
tained by FFT from H2 (two odd frames without object 
illumination), and typical location of the spatial filter 
(400 X 400 pixels). 
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Fig. 9. Reconstructed image from Hi i for a 400 x 400 
pixels (a) and 100 x 100 pixels spatial filter. 



without object illumination. The noise density is flat 
within the spatial filter region. Thus, the noise energy 
varies linearly with the spatial filter area. 

Since the signal increases quadratically with time, 
while the background noise increases linearly, the sig- 
nal to background ratio of the holographic detection in- 
creases linearly with time as expected for an ideal coher- 
ent detection. This mean that the equivalent noise is a 
fixed amount of energy that do not depends on time. On 
the other hand, since the signal is roughly flat with the 
area of the spatial filter, while the noise increases linearly 
with it, the signal to background ratio is proportional to 
the inverse of the spatial filter area. 

One must nevertheless avoid to reduce too much the 
filter area, since this area affects the image quality. This 
point is illustrated by Fig. [9]that compares reconstructed 
images of the USAF target made from Hi i with the 
slightly higher illumination level of Fig. HI The spa- 
tial filter is 400 x 400 pixels for Figl^a) and 100 x 100 
pixels for FiglQl^a). Since the area is larger on Fig. [H^a), 
the reconstruction involves more spatial frequency com- 
ponents yielding more noise, but also better resolution as 
seen by comparing Figs. [9] (a) and (b). 



C. Analysis of the noise 

As we have considered the case of an object (the USAF 
target) weakly illuminated, noise is not related to the 
signal of the object. It can therefore be measured in the 
absence of illumination of the object. This is what we 
have confirmed by the analysis of curves of Fig. [6l which 
shows that the intensity obtained with illumination 
in the black parts of the reconstructed image (bottom of 
the grooves associated with the image of the horizontal 
bars of the USAF target) is the same as obtained without 
object illumination. 

Without object illumination, the camera records only 
the attenuated signal from the laser (i.e. the LO refer- 
ence beam) , and the observed noise is equal to the sum of 
the technical noises with shot noise. The technical noises 
are mainly the electronic noise of the camera, the quanti- 
zation noise of the camera ADC (Analog to Digital Con- 
verter) and the laser intensity noise. Shot noise results 
from different stochastic processes (laser emission, ab- 
sorption by optical attenuators, photo conversion within 
the pixels of the camera) that yield an electronic signal 
corresponding for each frame and each pixel to an integer 
number of photo electron. This number is affected by a 
random Poisson noise (with a standard deviation equal 
to the root of the number of photo electrons). 

Let us consider the signal In{x^y) detected on frame 
n and pixel (x^y). Because of shot noise, In{x,y) is the 
sum of a statistical average component {I{x,y)) with a 
shot noise component o(x, y^ n), which is random for each 
pixel (x, y) and each frame n. Thus we have: 

In{x, y) = {I{x, y)) + On(x, y) (8) 

Since /n(^, y) is measured in photo electron units (e), the 
fluctuations Onix^y) obey Poisson statistic (or Gaussian 
since In{x^y) ^ 1) with : 

{on{x,y)) = (9) 
{\on{x,y)\') = {I{x,y)) 

Note here that Eq. ([9]) that defines On{x^y) is only valid 
if In{x,y) is expressed in photoelectron units (e). Fluc- 
tuations On{x^y) therefore depend on the gain of G the 
camera. The statistical average {In{x^y))^ and the shot 
noise On(x, y) terms of Eq. ([8]) are not accessible to mea- 
surement, and cannot be separated. Therefore, one can- 
not remove the shot noise term from a measured frame 
In{x,y). 

In the case of a camera, the photo electronic signal cor- 
responds to a relatively small number of photoelectrons 
e, and shot noise is relatively large. In our experiment for 
example, the local oscillator power is adjusted so that In 
corresponds to half saturation of the CCD camera. Since 
our camera is 12 bits. In is thus In ~ 2000 DC (Digital 
Count), i.e. In - 10"^ e (since G = 4.41 e / DC). The shot 
noise On{x^ y) ~ 100 e is thus much larger than electronic 
(30 e in our case) and quantization (4.4 e) noise. 

The relative importance of shot noise with respect to 
the technical noises depends heavily on the experimental 
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(a) 



Fig. 10. Principe of calculation of the Monte Carlo odd 
frames I'^{x^y) from Ii{x^y) recorded without object il- 
lumination. 



situation. In our case, the signal detected by the camera 
is dominated by the LO, and the signal is strong. The 
shot noise is important, and dominates the camera noise. 
In the case of a direct imaging (without interference with 
a LO) of an object weakly illuminated, the photoelectric 
signal is much weaker, and the electronic noise of the 
camera can become the dominant noise. 

To check if our holographic noise detection is or is 
not hmited by shot noise, we have calculated the odd 
frames by Monte Carlo simulation, and we have calcu- 
lated the corresponding noise. The main advantages of 
Monte Carlo simulation are that it is easy to implement 
and robust. Moreover, since the photons are emitted in a 
discrete and random way, Monte Carlo calculation yields 
intuitive description of our experiment, and of its noise. 

In our case, the detected signal is large: In{x^y) ^ 1. 
The statistical average {I{x^y)) is much larger than the 
fluctuation On{x,y). It can be thus replaced in Eq. (|9|) 
by the signal In{x,y) recorded on the same frame, or, 
since this signal is roughly the same for all frames, by 
the signal Ii{x^y) recorded without object illumination 
on frame 1. We have: 



{on{x,y)) = 
{\on{x,y)f) c:^ h{x,y) 



(10) 



One can calculate Onix^y) by Monte Carlo. Since all the 
holograms are built by making differences of frames (see 
EqlHand EqlTj), the contribution of the statistical average 
cancels in the holograms. As illustrated by Fig. [TOl one 
can thus replace {I{x^y)) by Ii{x^y) in Eql8]to calculate 
Monte Carlo odd frames I'^{x^y), and Monte Carlo odd 
holograms H'J 



^ni^^y) = h{x,y)^On{x,y) 



(11) 



The Monte Carlo odd holograms H'^{x^y) are similar to 
the measured odd holograms H'^{x^y). They are calcu- 
lated by replacing hni^.y) by I'^in+ii^.y) in Eq. [6l 

Figure [11] compares the reconstructed images from 4 
frames holograms measured without illumination of the 
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Fig. 11. Reconstructed images calculated from 4 frames 
recorded without object illumination i.e. from H'^ (a), 
and from 4 Monte Carlo frames i.e. from H'l (b). His- 
togram of the intensity \H\^ displayed in linear (black) 
and logarithmic scale (light grey): without object illumi- 
nation (c), and with Monte Carlo (d). The reconstruction 
is made with a 282 x 282 pixels spatial fllter. 



1e5 r- 



o1e4 




Fig. 12. Average intensity \H\^ as a function of the num- 
ber Ni, of frames: Ni, = 2 for H'^ or H!{, A^6 = 4 for 



H'^ or H'l,..., 



and Ni) = 32 for H^2 H^^. The frames 



32- 



are either recorded without object illumination H'^ or 
simulated by Monte carlo H'^. The images have been re- 
constructed with 2FFT. The size of the spatial fllter is 
400 X 400, 282 x 282, 200 x 200, 141 x 141 and 100 x 100. 



target (i.e. calculated from H^) in Fig lTTT a). with the 
Monte Carlo ones (i.e. calculated from H^) on Fig. 
fTTT b). The two images are very similar and show a Gaus- 
sian distribution for the holographic complex noise ampli- 
tude H{x,y,D), and an exponential distribution for the 
intensity \H{x,y, D)\'^ . Using ImageJ, we have plotted 
the histogram of intensity \H{x, y, D)\'^ for the measured 
(Fig lTTT c)) and the Monte Carlo reconstructed images 
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(Fig lTlT d)). The histograms plotted in logarithmic scale 
show that the measured and Monte Carlo distributions 
are exponential. This result is confirmed by the average 
intensity (measured: 2120; Monte Carlo: 2190) which is 
equal to the standard deviation (measured: 2121; Monte 
Carlo: 2196). We can also note that average intensities 
and standard deviations are the same for the measured 
and Monte Carlo data. The noise observed in experi- 
ments is thus equal to the shot noise that can be calcu- 
lated by Monte Carlo knowing the gain of G the camera. 

To confirm this results we have calculated the average 
intensities measured and simulated for holograms with 
different number of frames: A/'^ = 2 to 32, and for dif- 
ferent areas of the spatial filter: 100 x 100 to 400 x 400 
pixels. The results are presented in Fig(T2l Note that 
the measured points have already been plotted in Fig|71 
(a) (light gray curves). Figure [12] shows that, the mea- 
sured average intensities are extremely close to the Monte 
Carlo simulated ones for all numbers of frames, and all 
spatial filter areas. The noise observed in experiment is 
thus dominated by shot noise. 

It should be noted that thanks to the heterodyne gain, 
our setup is immune to electronic noises. Similar results 
are expected with any well-built scientific camera that 
uses low noise electronic in particular for the amplifying 
stage before analog to digital conversion. Industrial cam- 
eras that may exhibit electronic noise must be avoided. 
It should also be noted that the CCD exposure time (for 
1 frame) does not play any role, since it do not appear 
in the Monte Carlo simulation. Results depend only on 
the number of photo electrons. Thus, the exposure time 
plays no rule as far as the dark noise of the camera re- 
mains lower than the local oscillator signal (LO). 

The scaling laws observed for the experimental noise 
can be interpreted physically. The shot noise is a white 
noise, which is random and Gaussian in spatial and tem- 
poral frequencies. The holographic detection is coherent 
in space and time. As a result, on one hand (space), the 
detected noise is proportional to the number of spatial 
modes that are detected, i.e. to the number of pixels of 
the spatial filter. On the other hand (time), the noise is 
proportional to the product of the measurement time T 
by the detection bandwidth, which is here equal to 1/T. 
As a result, the equivalent noise does not depend on time, 
i.e. on the number of frames A/'^. 

One must notice that the holographic signal H{x^ 0), 
detected by the camera, is the interference of the optical 
field scattered by the object E with the reference beam 
local oscillator field Elo- It corresponds to E^qE. This 
signal is detected and integrated over time by the camera 
yielding a measurement in photoelectron units (e). This 
means that the reconstructed image, which corresponds 
to \H{x,y,D)\'^ (i.e. to |£^2oP x \E\'^ integrated over 
the measurement times), is measured in e^ units for each 
pixel. The vertical axis \H\'^ or of Fig. [5l Fig. [6l 
FigEand Fig IT2l corresponds to the product of the energy 
of the signal field \E\'^ in e units, by the LO field |£^loP 
in e units too. It is then possible to get an absolute 



calibration of the detected signal. 

Consider for example the image reconstructed from H^^ 
with a 282 x 282 spatial filter displayed on Fig (11] (a). In 
all our experiments, the frame signal In averaged over the 
pixels is In 6800 e. The LO signal is thus \Elo\^ =6800 
e for one frame, and 27200 e for 4 frames. Since is 
2120.53 e^ in that case (see Fig (11] (c) ), the equivalent 
noise \E\'^ is here 2120/27200 = 0.0779 e. This result is 
obtained for a 282 x 282 spatial filter. It corresponds to 
an equivalent noise of 0.0779 x (1024/282)^ = 1.02 1 
e without spatial filtering. This result is general and 
does not depends on the number of frames, nor on the 
spatial filter area. The equivalent noise corresponds thus 
1 photoelectron per pixel of the spatial filter whatever 
the number of frames or spatial filter area is^^^^^. 

It should be noted that the experimental and Monte 
Carlo simulated signals depends on the gain of the cam- 
era G with different factors. On one hand, the frame 
signal Inix^y)^ and the measured holograms H'^ are ex- 
pressed in photo electron units. They are proportional 
to G. On the other hand, the Monte Carlo simulation 
fluctuations On{x^y) are calculated from Eq. [TOj They 
depend on G through Iiix^y). The fluctuations On{x^y) 
varies thus like the square root of G. As result, the Monte 
Carlo holograms H'^{x^y) that are calculated from the 
Monte Carlo frames I'^{x^y) (defined by Eq. [11]), are 
proportional to y/G. 

Thus, the agreement of the measured noise with the 
Monte Carlo calculated shot noise depends on the ex- 
act value of the camera gain G. This gain G has been 
calibrated by the camera manufacturer (here G = 4.41 
e/DC) by illuminating the camera with a very clean white 
light source, which must be stable spatially and tempo- 
rally and by assuming that the observed fluctuations of 
the camera signal in time or in space result from shot 
noise^^. The holographic measurement we have done 
without object illumination is thus very close to what 
was done to calibrate the camera. Nevertheless, our light 
source is much less clean. Our laser may exhibit temporal 
fluctuations, and, because of stray reflections, the power 
of the LO beam is not constant on the camera detect- 
ing area, as shown on image of frame Ii at left in FiglTOl 
The good agreement between the measured noise and the 
noise calculated by Monte Carlo from the gain G of the 
camera, shows that the fluctuation we got with a laser 
are the same than the fluctuation observed by the cam- 
era manufacturer with a white light source yielding the 
same average camera signal. This means that the defects 
of our experimental setup (temporal fluctuations of the 
laser intensity and spatial variations of the LO signal) 
do not seem to bring excess of noise with respect to shot 
noise. 

We think we are insensitive to defaults of our setup 
because we do both heterodyne and off-axis holography. 
The holographic signal is thus modulated both in time 
(heterodyne) and space (off-axis). The unwanted low fre- 
quency time components are filtered off by building holo- 
grams with many frames (see Eq. [6] and Eq. [7]), while 
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the low frequency space components are filtered off by 
spatial filtering. 

4. Conclusion 

By recording hologram of a US AF target under weak illu- 
mination with our off axis heterodyne holography setup, 
we have found that the signal and noise scales as ex- 
pected with a shot noise limited experiment. The signal 
for the reconstructed image intensity \H\'^ is proportional 
to the square of the time, while the noise, which is equal 
to the expected shot noise^^, is proportional to the time, 
and to the number of modes involved in the holographic 
reconstruction (area of the spatial filter in pixels units). 
These results validate previous weak holographic signal 
experiments^^"— 
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